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ABSTRACT 

The first mechanism invoked to explain the existence of the thick disk in the Milky Way Galaxy, were 
the spiral arms. Up-to-date work summon several other possibilities that together seem to better 
explain this component of our Galaxy. All these processes must affect differently in distinct types 
of galaxies, but the contribution of each one has not been straightforward to quantify. In this work, 
we present a first comprehensive study of the effect of the spiral arms in the formation of thick 
disks, as going from early to late type disk galaxies, in an attempt to characterize and quantify this 
specific mechanism in galactic potentials. To this purpose, we perform numerical simulations of test 
particles in a three-dimensional spiral galaxy potential of normal spiral galaxies (from early to late 
types). By varying the parameters of the spiral arms we found that the vertical heating of the stellar 
disk becomes very important in some cases, and strongly depends on the galaxy morphology, pitch 
angle, arms mass and its pattern speed. The later the galaxy type, the larger is the effect on the 
disk heating. This study shows that the physical mechanism causing the vertical heating is different 
from simple resonant excitation. The spiral pattern induce chaotic behavior not linked necessarily to 
resonances but to direct scattering of disk stars, which leads to an increase of the velocity dispersion. 

We applied this study to the specific example of the Milky Way Galaxy, for which we have also added 
an experiment that includes the Galactic bar. From this study we deduce that the effect of spiral 
arms of a Milky-Way-likc potential, on the dynamical vertical heating of the disk is negligible, unlike 
later galactic potentials for disks. 

Keywords: galaxies: evolution - galaxies: kinematics and dynamics — galaxies: spiral — galaxies: 
structure 


1. INTRODUCTION 

Simulations of galaxy formation are coming to point 
where detailed processes of galaxies have never been 
explored before in detail, such as random and rota¬ 
tional velocities can be better studied and understood 
dScannapieco ct aDl20lL l . Details on disk potentials can 
be probed and compared with observations, and we are 
now able to shed some light on evolution of galaxies start¬ 
ing now on small scale stellar motions. 

Dynamical heating of the Milky Way disk has now been 
known for over 60 years, mainly through observations in 
the solar vicinity. From those observations we learnt that 
stellar random motions corr elate nicely with their ages 
know n as the age-cr relation (lWielenlll977t iBinnev et al.1 
I2000H . In particular, in the case of the Milky Way disk, 
it is known that the radial velocity dispersion is twice as 
much as the vertical dispersion and that the radial scale 
length of the thi ck disk is muc h shorter than that of the 
thin disk (|Bensbv et al. l[20TT1i . 

Recent studies also show that many, if not all, edge- 
on spira l ga laxies appear to h ost dual disk systems 

dGerssen fc Shapiro Griffmll2012llva,n der Kruit fe Searlel 

a younger, dynamically colder and thinner com¬ 
ponent: the thin disk and at least one older compo¬ 
nent, (mainly) stellar, dy namica lly hotter and thi cker 
component: the thick disk (lYoach im fc Dalcanton 1120061 

lYoachim fc Dalcantonll2~008HComer6n et al.1 20111b 

It is still not straightforward to dilucidate the 


mechanism responsible for the vertical heating of the 
disk from observations, specially since more likely it 
might be rather a combination of several possibil¬ 
ities. Among the mechanisms proposed there are 
some external to the disk such as hits by satel- 


19971; IVelazauez fc White: 19991 Benson et al.l 2004; 

Font et al.l 120011 Quinn et al. 

19931 iVillalobos & Helmi 

20081: Di Matteo et al.l 2011 

); scattering by dark 

halo obiects or globular clu 
20091: Hanninen & Flvnnj |20C 
ternal origin, such as, dyi 
rect encounters with giant m 

sters dVande Putte et al. 
"2|). And some of in- 

lamical heating by di- 

olecular clouds (iGarlbera 

1987 

IVillumsenl 19831 Lacev & Ostriker 19851 Lacev 

1984 

: iSoitzer & Schwarzrchildl 19511 Inouc fc Saitoh 


2014); heating by encounters with th e potential 


produ ced by long-lasting spiral arms dFaure et al.l 
12011 or irregular and transient spiral structure 

(iMinchev fc QuilIeiill2006t lFuchsl[20(Htl,Ienkins fc Binnevl 

ll990HBarbanis fc Woltiedl l967fl: pertur bations from stel¬ 

lar bars ( Saha. Tseng fc Taaml I2010H : dissolution of 
young stellar clusters (|Kroupal l2002fl : or during an in¬ 
tense star formation phase in a period of i ntense accre¬ 
tion v ery early in the history of the Galaxy dSnaith et al.l 
(201411 . 


Because of the nature of these theories, the effects are 
dependent on galaxy morphology, particularly the intrin¬ 
sic mechanisms such as bars and spiral arms. Therefore, 
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to deeply understand secular evolution of disk galaxies, 
it is critical to study dynamical heating in a good sample 
of different disk galaxy types. Finally, from the point of 
view of observations, the radial heating agent seems to 
vary exactly as expected if the agent were the spiral arms, 
which provides a good chance that the spiral structure 
has at least an important rol e as a heating mechanism 
i n th e plane of galactic disks dGerssen fe Shapiro Griffinl 
I2012T) . The importance and influence of spiral arms, 
and even their very same nature is still under debate, 
there is no straightforward observational prove yet of 
their effect on stars, however it is nowadays consid¬ 
ered to_have a_kev role on large-sc ale galactic dynam¬ 
ics dSellwoodl [201311 for a review, (IRoskar et ali l2012t 


Minchev et al. 


2012t iLepine et all 120111 : iQuillen et ahl 


20111 : lAntoia et al.l 12009th One plausible observational 


space, known as “moving groups” dProctoilll869l: lEggen 

1959, 

1977L 199(1 1996allb 

; 1 Wilson et al. 19231; Roman 

1949; 

Soderblom & Mavoi 

19931 IMaiewskil" 1994: 1996 


rect though, eviden c e of t he effect of the spiral arms 
fChe reui et all [19981 [l999i: [Pe hnenl fl998t iFamaev et alJ 

120051120081: lAntoia et al .11201 2HPomneia et al.ll201llf . Of 


course, the spiral arms are not the only nor the prefered 
mechanism to explain moving groups, the galactic bar is 
other possibility. 

In this work we focus on the very firs t proposal to 
explain the vertical heating dWiclenl 11977 b disregarded 
at some point in the history because of their negligi¬ 
ble effect on the Milky Way disk: the dynamical heat¬ 
ing by effect of the spiral arms. We attempt to isolate 
and quantify the contribution of the spiral arms to the 
disk heating of galaxies. We performed numerical sim¬ 
ulations of test particles in a thre e-dimensional galacti c 
potential that models spiral arms ([Pichardo et alil2003 ), 
adjust ed to simulate spiral galaxies, f rom early to late 
types (jPcrez-Villegas et~aLll2012t 1201314 . We produced a 
set-up with relaxed initial conditions for a stellar disk. 
Finally we calculated the effect on the vertical heating of 
the stellar disk produced by the nonaxisymmetric large 
scale structures. We have included a preliminary study 
on a self-gravitating potenti al, known as PERLAS model 
(IPichardo et al.l[200311200411 for the Milky Way Galaxy, 
that includes the spiral arms and the galactic bar. 

This paper is organized as follows. The galactic mod¬ 
els, initial conditions and methodology are described in 
Section [2] The role of each one of the parameters of 
the model is studied with detail in Section [3J where we 
present calculations of dispersion velocity, the velocity 
ellipsoid, time evolution for spiral galaxies from early to 
late types, and an application to the Milky Way Galaxy. 
Finally, in Section 01 we present our conclusions. 


2. METHODOLOGY AND NUMERICAL IMPLEMENTATION 

The effect of the spiral arms over the stellar disk 
has b een studied profusely in either N-body s i mula¬ 
tions (ISe llwood fe Binnevll2002t IRoskar et al.ll2012l 120131 
iKawata et al.l 1201414 or with spiral patterns treated as 
perturbations to the axisymmetr ic background and mod¬ 
ellin g it as a density wave El Simone et ahl 12004 
iMinchev fe Quillenl 120061 : iFaure et al.ll2014lj~ 

In N-body simulations, although self-consistent, it 
is not plausible to adjust a given specific galaxy, or 


to isolate the effect for example, of the arms or es¬ 
tablish in detail the role that each one of the pa¬ 
rameters that characterize the spiral pattern play over 
the vertical heating of the disk. On the other hand, 
when treated as steady spirals exerted as a perturba¬ 
tion to the a xisymmetric poten tial, according to the 
hypothesis of iLin fe Shu I (1196914 . the heating is mini¬ 
mum a nd linked only to the resonant regions of the 
spira ls dLvnden-Bell fe Kalnaisl 119721 : iDe Simone et al.l 
1200411. being more efficient in the radial directio n 
(ISellwood fc Garlberg lfl984 : Minchev fc Quillenl[2006 1). 

Here we use spiral arms that although steady, are very 
different in nature to the typically and widely employed 
in literature. The gravitational potential due to the spi¬ 
ral pattern is not a simple perturbation but is rather 
based on a mass density distribution. With this model 
of spiral arms our studies on the vertical heating of the 
disk contrast in general with the density wave approach 
(except of course, for really small spiral arm masses and 
pitch angles). 

For the orbital study we employed then three- 
dimensional galactic potentials to model normal spiral 
galaxies (Sa, Sb and Sc). The motion equations are 
solved in the non-inertial reference system of the spiral 
arms and in Cartesian coordinates (x ', y', z'). The orbits 
are integrated for 5 Gyr with a Bulirsh-Stoer algorithm 
(Press et al. 1992), with a conservation of the Jacobi 
constant approximately up to 10 -12 . The disk heating is 
computed through the measure of the velocity dispersion 
at different times. 

2.1. Models for Normal Spiral Galaxies 

The models include an axisymmetric component 
(bulge, massive halo and disk), as the background po¬ 
tential, formed by a Miyamo to fc Nagai (1975) dis k 
and bulge, and a massive halo (lAllcn fc Santillanlll991l f. 
The parameters used to model normal spiral galaxies 
(Sa, Sb and Sc) are prese nted in Table Q] (compiled by 
iPerez-Villegas et all (1201314 4. 

Superposed to the axisymmetric components, for the 
spiral arms potential, we employed a bisymmetric self- 
gravitating three-dimensional potential, based on a den¬ 
sity d istribution, called PERLAS model (IPichardo et al.l 
1200314 . This potential consists of individual inhomo¬ 
geneous oblate spheroids s u perposed along a logarith¬ 
mic spiral locus ([Roberts. Huntley fc van Albada 111979( 4. 
Each spheroid has a similar mass distribution, i.e., sur¬ 
faces of equal density are concentric spheroids of con¬ 
stant semiaxis ratio. The model considers a linear fall 
in density within each spheroid. The minor and major 
semiaxes of each oblate spheroid are 0.5 and 1.0 kpc, re¬ 
spectively (this gives a width of the spiral arms of 2 kpc 
and height of 0.5 kpc from the disk plane) and the sepa¬ 
ration among the spheroid centers along the spiral locus 
is 0.5 kpc. The superposition of the spheroids begins 
and ends, in the ILR and CR, respectively. The density 
falls exponentially along the spiral arm, where the radial 
scale length of the galactic disk is used (depending on 
morphological type, see Table 1). 

The mass assigned to build the spiral pattern is sub¬ 
tracted from the disk mass to keep the given model in¬ 
variable in mass. PERLAS is a more realistic potential 
since it is based on a density distribution and considers 
the force exerted by the whole spiral structure, obtain- 
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Table 1 

Parameters of the Galactic Models 


Parameter 


Value 


Reference 


Sa 

Sb 

Sc 


Axisymmetric Components 

M b / Mo 

0.9 

0.4 

0.2 

1,2 

Mfl/Mfl 1 

0.07 

0.09 

0.1 

2,3 

Rot. Velocity (kms -1 ) 

320 

250 

170 

4 

(y D /v Rot )2 Rh 

0.51 

0.65 

0.70 


M d (10 n M Q ) 

1.28 

1.21 

0.51 

3 

M s (10 11 Mq) 

1.16 

0.44 

0.10 

Mb/Mjj based 

Ms (10 n M Q ) 

16.4 

12.5 

4.8 

Mb/Mb based 

Disk scale-length (kpc) 

7 

5 

3 

1,3 

bi 2 (kpc) 

2.5 

1.7 

1.0 


a 2 2 (kpc) 

7.0 

5.0 

5.3178 


b 2 2 (kpc) 

1.5 

1.0 

0.25 


a 3 2 (kpc) 

18.0 

16.0 

12.0 


Spiral Arms 

locus 


Logarithmic 

5,9,10 

arms number 


2 


6 

pitch angle (°) 

8-40 

9-45 

10-60 

4,7 

M sp /Ms 


1-5% 


9 

scale-length (kpc) 

7 

5 

3 

disk based 

(km s -1 kpc -1 ) 

-30 

-25 

-20 

5,8 

ILR position (kpc) 

3.0 

2.29 

2.03 


CR position (kpc) 

10.6 

11.14 

8.63 


inner limit (kpc) 

3.0 

2.29 

2.03 

~1LR position based 

outer limit (kpc) 

10.6 

11.14 

8.63 

'-'-'CR position based 


References. — (1) Weinzirl et al. 2009; (2) Block et al. 2002; 

(3) Pizagno et al. 2005; (4) Brosche 1971; Ma et al. 2000; Sofue & 
Rubin 2001; (5) Grosb0l & Patsis 1998; (6) Drimmel et al. 2000; Grosb0l 
et al. 2002; Elmegreen & Elmegreen 2014; (7) Kennicutt 1981; (8) Patsis 
et al. 1991; Grosb0l & Dottori 2009; Egusa et al. 2009; Fathi et al. 2009; 
Gerhard 2011; (9) Pichardo et al. 2003; (10) Seigar & James 1998; Seigar 
et al. 2006; 

1 Up to 100 kpc halo radius. 

2 bi, a 2 , b 2 , and a 3 are scale lengths. 

ing a more detailed shape for the gravitational poten¬ 
tial, unlike a two-dimensional local arm such as the tight 
winding approximation (TWA) represented for a simple 
cosine function. 

Spiral arms nature is a matter of discussion nowadays, 
particularly their long-lasting or transient nature. We 
have performed experiments with constant, transient, 
gradual and sudden presence of the spiral arms. Al¬ 
though the growth rate is an unknown parameter in 
galaxies, we have considered different cases to test. On 
one hand we produce a set of experiments where the to¬ 
tal mass of the spiral arms is introduced at once (t = 0 
Gyr). The second set of experiments, inserts the spiral 
arms mass linearly in a timelapse of 1 Gyr. And a third 
set of experiments for which the spiral arms are simu¬ 
lated as transient, they vanish and grow with a given 
periodicity. 


2.2. Initial Conditions and Equilibrium of a Stellar Disk 

The initial conditions set-up follows the Miyamoto- 
Nagai density profile we are imposing. This to avoid 
transient effects induced by differences between the ini¬ 
tial particle distribution and the imposed disk potential. 
In this manner, the initial condition for the stellar disk 
is given by 


_ b\M d a 2 R 2 + (a 2 + 3a/ z 2 + b\)(a 2 + \Jz 2 + 6|) 2 

PMN ~ ~4^ (R2 + ( ffl2 + ^ z 2 + &2)2)5/2( z 2 + ^ 2 ) 3/2 

(!) 

where M d is the mass of the galaxy disk, a 2 and b 2 are 
the radial and vertical scale-length, respectively. This 
three parameters span a range of values in our simula¬ 
tions in order to capture different galaxy morphologies 
and kind of spiral arms. 

To distribute the particles according to the Miyamoto- 
Nagai density law we solved equation © with a root 
finder method. This is done by expressing the den¬ 
sity p(R,z) in terms of the ratios p(R,0)/p(0,0) and 
p{R , z)/p(R, 0), which provides us an equation for R and 
z in terms of the density. The value of these ratios ranges 
from 0 to 1, therefore we can explore all the possible val¬ 
ues of the density with a random function and solve it 
for R and z. 

To assign veloci ties to the particl es we follow the strat¬ 
egy proposed by iHernauisH •: 11)93 : where velocities are 
distributed by an approximation using moments of the 
collisionless Boltzmann equation plus the epicycle ap¬ 
proach. 

Thus we proceeded as follow: once the density pro¬ 
file has been established, it is necessary to obtain the 
rotation velocity. This can be derived from <f>, the grav- 
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itational potential of the model 

/ 1 <9<F\ 1/2 

< 2 > 

and v c = RQ C (R), so the circular velocity is given by 

M R )=(R M ) • (3) 

Once known f2 c at any radii we obtain 

/ wf } 2 \ 1 / 2 

*=( 4!! ‘ + J ^) ' (4) 

known as the epicyclic frequency, necessary to calcu¬ 
late the velocity dispersion at R and to correct for the 
asymmetric drift. 

To achieve the requirement for the stellar disk to be 
in equilibrium, it is necessary to introduce a given dis¬ 
persion in the velocity as a function of R. The velocity 
dispersions in the three polar coordinates are 


„* = 3.368 E(B)Q 

(5) 

K 

1 OrR 

(6) 

^ ~ 2 n c 

v z = ^/nGTjb 2 , 

(7) 


where k is the epicyclic frequency, E(I?) the surface den¬ 
sity, &2 is the vertical scale-lenght of the disk, and Q 
is the known Toomre parameter. According to iToomrd 
(1964!), local stability requires Q > 1, we choose Q = 1.1 
and found this value to be sufficient for the three galaxy 
types. In this way the velocity dispersion depends on the 
mass of the components that form each galaxy. 

The asymmetric drift correction is defined as 

M 2 = v l ^ 

is a correction that has to be implemented in the set up 
for the initial conditions given the fact that stellar orbits 
are not in general in circular orbits, instead, orbits follow 
epicycles around a guiding point at the position of the 
circular orbit and these epicycles are characterized by the 
epicyclic frequency k. 

Finally the particles are distributed in the velocity 
space as 

= (v^) ± x<70 

VR = ( Vr ) ± XCTr (9) 

v z = (v z ) ± xa z 

where a; is a random number between 0 and 1, (v#) is 
given by Eq. © and the average radial and vertical 
velocities are taken as (vr) = (v~) = 0. 


2.3. Dispersion analysis 

The disk heating is often referred as the increase in the 
velocity dispersion over the lifetime of a star. Any disk 
thickening is then related to an increase in the vertical 
velocity dispersion of the disk stars. In this paper we 
analyze the spiral arms effects on the stellar disk, based 
on the study of the vertical velocity dispersion and its 
dependence with the parameters that characterize the 
spiral pattern. 

The vertical velocity dispersion tr 2 , is then calculated 
in the simulations by dividing the plane z = 0 into 1 kpc 
bins and computing, as usual, the squared root of the av¬ 
eraged squared vertical velocity for all the particles that 
fall into a given bin. This provide us the vertical veloc¬ 
ity dispersion as a function of R. In order to establish 
the contribution of the spiral arms to the disk thickness 
we also compute a z as a function of time by measuring 
the velocity dispersion at a fixed radius R for every time 
code unit (in this case, every 100 Myr across 5 Gyr in the 
simulation). 

2.4. Control Simulations: Testing the Initial Conditions 
Equilibrium 

As described in section c m the first goal is to build 
an initial stellar disk in equilibrium to be sure that any 
change seen in radial and vertical velocity dispersion is 
strictly due to the interaction of the spiral arms with the 
stellar disk and not originated by a spurious non-relaxed 
initial condition set-up. 

For this test a control simulation was produced with 
only the axisymmetric components for the potential 
model. Stars are run by 5 Gyr, for all galaxy types. Fig¬ 
ure H] shows a, as a function of R at different stages in 
the temporal evolution for an Sa, Sb, and Sc galaxy. 

From the figure it is clear that the vertical velocity 
dispersion do not evolve or deviate from the initial dis¬ 
persion, as expected for a disk in equilibrium with the 
axisymmetric potential. 

3. RESULTS 

We present in this section a set of controlled experi¬ 
ments to study the dynamical heating of disks on spiral 
galaxies considering the spiral arms as the driver. The 
general purpose is to shed some light on the relative im¬ 
portance of these large-scale structures to other sources 
of dynamical heating in different morphological types. 
The experiments include studies of different structural 
and dynamical parameters of the spiral arms such as, 
different pitch angles, total spiral arms masses, angular 
speeds, transient, and one final case modeling the Milky 
Way Galaxy (with preliminary results that will be better 
developed in a future work). 

3.1. Dependence of the disk heating induced by spiral 
arms with galaxy morphology 

Spiral galaxies present a wide variety on morphological 
types, from massive bulge-dominated galaxies to practi¬ 
cally bulgeless disks, spanning a wide range of values for 
the parameters that characterize different galaxy types. 

Figure [2] shows the velocity dispersion, a z , as a func¬ 
tion of the galactocentric radius, R , at different times in 
the simulation for our galactic models: Sa, Sb and Sc 
(introduced in Section [2]). For these three simulations 
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2 time(Gyr) 3 





Figure 1. Test for the initial conditions equilibrium: velocity 
dispersion of the stellar disk in the axisymmetric potential only, as 
a function of R , along 5 Gyr time evolution for each galaxy type. 

the mass of the spiral arms, M arms , is 5% of the total 
disk mass with a Pitch angle of 40° for the Sa galaxy, 45° 
for the Sb galaxy, and 40° for the Sc galaxy. For these 
experiments we have employed the largest spiral arms 
masses and pitch angles, for plausible (non-fully chaotic) 
galactic models, to identify clearly spiral arms effects, if 
any. The three plots show a distinct increase in the ver¬ 
tical velocity dispersion caused by the spiral arms. Also 
from figure[2]it is clear that the change in a z with respect 
to the initial dispersion is smaller for the Sa galaxy and 
grows with the morphological type, being much larger 
for the Sc galaxy model. 


2 time(Gyr) 3 



Figure 2. As in Figure [l] but including now the spiral arms 
contribution to the potential. The vertical velocity dispersion is 
plotted as a function of R for 5 Gyr time evolution. 

The dependence of the effect of the spiral arms with 
the morphology is such that for an Sc galaxy the effect 
is evident in the spatial distribution of the stellar disk 
particles. Figure [3] shows the x — z projection of the 
stellar distribution plotted at t = 0, t = 2.5 Gyr, and 
t = 5 Gyr. 

Also, a thickening of the disk is discernible during the 
orbital evolution when comparing with the initial distri¬ 
bution. 

The thickness can be quantified by computing the root 
mean square of the coordinate z, i.e. z rms = y/< z 2 >. 




















6 


Martinez-Medina et al. 2015 



x(kpc) 


Figure 3. Edge on disks show the evolution of stellar orbits at three different times (0, 2.5 and 5 Gyr). The model corresponds to an Sc 
galaxy. The thickening of the disk due to the spiral arms presence is clear. 



Figure 4. Disk thickness z r ms as a function of R for t = 0, 2.5 Gyr 
and 5 Gyr. 

Figure [4] shows the thickness as a function of R for the 
stellar disks in Figure [3] 

From these first set of experiments, we conclude that 
the sharpest effect on the velocity dispersion, is present 
on the latest morphological types. It is worth mention¬ 
ing here that, although we are separating the models in 
Sa, Sb and Sc galaxies, with strong gaps in between the 
different models in the initial scale height, the results on 
the disk heating driven by spiral arms presented here, are 
more general, i.e. the observed heating results significant 
in thinner disks, which, as a consequence, has implica¬ 
tions on galaxy types, in this case, particularly on later 
types. 

For the earliest type (Sa), the isolated effect of spi¬ 
ral arms corresponds to a maximum increment of 7% 
percent of the initial velocity dispersion, this consider¬ 
ing the most massive and the largest pitch angles pos¬ 
sible to produce plausible galactic models. Likewise, for 


the intermediate galactic type (Sb), the isolated effect 
of spiral arms corresponds to a maximum of increment 
20% percent of the initial velocity dispersion, again con¬ 
sidering the most massive and the largest pitch angles 
possible to produce plausible galactic models. Finally, 
for the latest type (Sc), the isolated effect of spiral arms 
corresponds to a maximum of 62% percent of the ini¬ 
tial velocity dispersion for this example, where we have 
not used the maximum plausible parameters for the spi¬ 
ral arms. When the maximum pitch angle and mass is 
employed, the percentage goes up to almost 90% of the 
initial velocity dispersion. 

Since a visible effect in the vertical dispersion when 
spiral arms are included is considerably larger for the 
latest galactic types, compared with early and interme¬ 
diate types, in the next sections we concentrate on a more 
detailed study of the disk thickening focusing only in the 
late type galaxies. 

3.1.1. The pitch angle effect 

The pitch angle is one of the most influential structural 
parameters that characterize spiral patterns. In this sec¬ 
tion a range of values is explored in order to quantify 
the dependence of the increment in the vertical velocity 
dispersion with pitch angle in the most affected galactic 
models, that are the latest types. 

Figure [5] shows three plots of a z vs R , where the mass 
of the spiral arms is set to a constant for each plot and the 
pitch angle vary according to Table 1. Each plot has the 
initial dispersion curve a z (R,t = 0) and the dispersion 
after a 5 Gyr evolution a z (R, t = 5 Gyr) for each pitch 
angle value. 

From Figure 0 it is clear that regardless the mass of 
the spiral arms, the vertical velocity dispersion increases 
notably with the pitch angle. The less massive the spiral 
arms the smaller is their effect in general, as expected. 
Indeed, for spiral arms masses smaller that ~ 1% of the 
disk, the contribution of spiral arms to the dynamical 
heating becomes negligible. 
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Pitch Angle 



R(kpc) 


Figure 5. Final velocity dispersion after a 5 Gyr evolution as a 
function of radius for spiral arm masses of 1%, 3%, and 5% of the 
total disk mass. 


velocity dispersion. We find that this occurs at approxi- 
mately R m ax = 1-5 kpc for this model. We then keep R 
constant at that value and plot A cr z there as a function 
of the pitch angle. This will show us the tendency seen 
in the previous section. Now repeating this for the three 
spiral arms masses employed will show how A a z scales 
with this second parameter. 

Figure [6] shows Aa z as a function of the pitch angle for 
three different masses of the spiral arms. Here it is clear 
that both parameters, the pitch angle and the mass of 
the spiral arms, affect considerably the disk thickening 
effect driven by the spiral structure. 



Figure 6. Difference in the vertical velocity dispersion Acr z (at 
t=0 and t=5Gyr) vs the pitch angle, plotted for the three masses 
of the spiral arms, as a scaling factor. 

We produce with these results an empirical functional 
relation between Act, and the Pitch Angle. The fit of 
data plotted in figure [6] is made by noting that A a z , 
increases slowly at small angles, then the slope of the 
curve grows with the angle and flattens for the largest 
values of the pitch angle. This behavior could be inter¬ 
preted as a saturation effect of A a z after a certain time, 
this time being shorter for higher angles and masses. 
The s aturation effect is also seen in a a - time re¬ 
lation dSeabroke fc Gilmorell2007t iSoubiran et alJ [20081 
iCalberg et al.lll985fl . where the dispersion remains con¬ 
stant after ~ 5 Gyr. Based on this observation, a nice 
fit to the results would be a Boltzmann sigmoidal func¬ 
tion which is characterized by displaying a progression 
from small beginnings that accelerates and approaches 
a climax over the independent variable. The Boltzmann 
sigmoidal function, for this particular case, is defined by 


3.1.2. The spiral arm mass effect 

As it is shown in the previous section, the effect of the 
pitch angle can be significant to the disk thickening. It 
also results intuitively clear that this also scales with the 
mass of the spiral arms. To address this point, we pro¬ 
duced several experiments condensed in Figure [Gj that 
show the disk thickening dependence on the spiral arms 
mass and the pitch angle together. For this purpose, 
first we identified the radius Rmax at which occurs the 
maximum difference between the final and initial vertical 


Aa~ 


A 2 + 


{A\ — A2) 

1 + exp((x — Xo )/d) 


where: 

x = Pitch Angle (°) 
xo = center (°) 
d = width (°) 

A\ = initial Act- value (kms -1 ) 
A 2 = final Ao z value (kms -1 ) 


(10) 
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The center Xq is the pitch angle at which A a z is 
halfway between A\ and A 2 . The width d is related with 
the steepness of the curve, with a larger value denoting 
a shallow curve. 

Figure [7] shows the fit of the sigmoidal function (Eq. 
EO) to the data. We see that this function reproduces 
well the behavior of A cr z with the pitch angle for the 
three masses of spiral arms used in our simulations. 
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Figure 7. Fit of the Acr z - pitch angle relation with the Boltz¬ 
mann sigmoidal function for the three masses of the spiral arms. 

Taking Eq. m as the functional form for the A<r~ - 
pitch angle relation, in Table [21 we summarize the pa¬ 
rameters that describe the fits plotted in figure [7] for the 
three spiral arms masses. 


Gradually Increasing the Mass of the Spiral Arms 

As explained before, the spiral structure in our galaxy 
potential models is imposed and fully introduced since 
the starting point of the simulation. However, one might 
wonder if spurious effects (such as a drastic dispersion 
increase in the early stages of the simulation) are intro¬ 
duced with this method; also, in real galaxies the birth 
and death of spiral arms is most probable not a sudden 
process. 

In order to explore these two scenarios, we prepared 
a set of representative simulations modifying the model 
in such a way that the arms are allowed to have a grow¬ 
ing period until the total assigned mass is reached. For 
this set of simulations we started with an Sc galaxy 
with different combinations of values for the spiral arms 
mass and the pitch angle. We chose two radial positions 
and measured in there the evolution of a z with time as 
showed in Figure[8] In this figure the temporal evolution 
of cr z for the sudden full mass arms is compared with a 
model with linear (1 Gyr period) growing mass. It is clear 
that there is a shift in the dispersion achieved at the end 
of the simulation, but most important, for both kinds of 
imposed spiral arms the velocity dispersion increases at 
roughly the same rate. 

With this simple experiment we are not pretending to 
capture the great variety of processes that lead to the 
formation of a real spiral pattern, nor their timelapse of 



Figure 8. a - t relation that shows, regardless the radial position, 
a continuous heating for both kind of imposed spiral arms: those 
that grow with time (blue lines) and those totally formed since the 
beginning (black lines). 
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Table 2 

Act, - Pitch Angle Fit Details 


Marms {%Mdisk) 

xo (°) 

d(°) 

Ai (km s 1 ) 

A 2 (km s 1 ) 

i% 

64.73 

24.05 

-0.33 

5.32 

3% 

54.73 

14.95 

-0.057 

24.84 

5% 

35.68 

6.82 

1.69 

28.44 


Note. — Parameters that combined with equation 1101 describe 
the A(7 z - pitch angle relation. Center xo, width d, initial Acr z value 
Ai, and final A a z value A 2 . For the three spiral arms masses (in 
percentage of the disk mass Mdisk)- 


full action, but this is still likely a better approximation 
than a sudden action of spiral arms. With this excer- 
cise we notice however that, for any radial position, the 
vertical dispersion is only shifted when using spiral arms 
that have a growing period. This do not suppose impor¬ 
tant differences with our previous simulations or results 
derived from them. 

3.1.3. Varying the Angular Velocity of the Spiral Pattern 

The vertical heating of the stellar disk produced by the 
spiral arms depends on the parameters of the nonaxisym- 
metric structure, particularly for the late galaxy types. 
By varying the pitch angle and the mass of the arms we 
have found a correlation between this parameters and 
the thickness of the stellar disk. 

In this section, we study now the pattern speed effect 
on the disk vertical heating. For this purpose, we ran 
a set of simulations with different values of the spiral 
pattern for the latest types of galaxies: Sb and Sc, that 
are more clearly affected for the structural parameters 
of the arms. I11 Figure [9] we show first the evolution of 
the vertical velocity dispersion <j z over time up to 5 Gyr 
for different values of the pattern speed, D, for an Sc 
galaxy. There seems to be a clear relation between the 
angular velocity and the dynamical heating. The velocity 
dispersion increases for slower rotating patterns. 

Figures [TO] and [Tl] summarize the results of the set of 
simulations when varying the pattern speed for Sb and Sc 
galaxies. The final stage of <j z is shown across the entire 
disk for each one of the pattern speeds used. As expected, 
the final vertical velocity dispersion and disk thickness 
are larger for slower rotating arms, independently of the 
galaxy type. Smaller values of D allow the spiral arms 
to heat more efficiently the stellar disk, likely due to the 
minor relative angular velocity between the arms and the 
stars, this allows the stars to interact with the potential 
of the arms for longer periods of time. 

Other studies link the heating of the stellar disk 
and the role of the patt ern speed on ly to the reso¬ 
nant regions of the spi rals (|Lvnden-Bell k. Kalnaislfl972t 
iDe Simone et al.ll2004lh in contrast, in this study we find 
that the heating, radial and vertical, occurs along the 
entire length of the spiral arms. 

On the other hand N-body simulations had shown to 
develop transien t spira l structure t h at spans a range of 
patte rn speeds (|Sellwood fe Binnevl 120021 IRoskar et al.1 
l2012f) , but as the arms are transient and is not possible to 
isolate its effect in N-body simulations, it results difficult 
to establish a dependence of the heating on the pattern 
speed. 


£2(km ^ kpc" 1 ) 



Figure 10. Final velocity dispersion after a 5 Gyr evolution for 
different spiral pattern angular speeds in an Sb galaxy. 


£}(km i'J kpc' 1 ) 



Figure 11. Final velocity dispersion after a 5 Gyr evolution for 
different spiral pattern angular speeds in an Sc galaxy. 

The experiments and the model shown here allows us 
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Figure 9. Time evolution of the vertical velocity dispersion for three different spiral pattern angular speeds in an Sc galaxy. This 
comparison shows that slower rotating spiral arms heat more efficiently the stellar disk. 


to establish a relation between the vertical heating of the 
stellar disk and the pattern speed. 

3.1.4. <t z - Time Relation 

It is already known that the age and velocity dis¬ 
persion of stars are correlated. This has been estab¬ 
lished from observations in the solar neighbourhood 
as we l l as from numerica l simul ations dHolmberg et all 
mm IRoskar et all l20Tl 120141 Q). The a - t rela¬ 
tion shows a smooth, general increase of the veloc¬ 
ity dispersion with time and is best parametrized by 
a power law with exponents ran ging between 0.2 - 0.5 
(jGerssen fc Shapiro Griffinll2012l l. 

We explore the er - t relation in our simulations to find 
out if the velocity dispersion in the stellar disk due to the 
spiral arms fits with a power law t a and more important, 
to establish a range of values for a. 

To measure the time evolution of a z in our simulations 
first we locate the radius R m ax at which occurs the max¬ 
imum increase in the vertical velocity dispersion. This 
radius is Rmax = 1.5 kpc and is the same for all the sim¬ 
ulations with the Sc galaxy, independently of the pitch 
angle or the mass of the spiral arms. 
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log(a z ) dependency. The most of the points fall on a 
straight line, this reveals a behavior of the form cr z oc t a . 

Figure [13] shows the evolution of a z with time at 
R = Rmax , M arms /M d i sk = 0.05 and at different pitch 
angles: 20°, 30°, 40° and 50°. The black line in each plot 
is the best fit of the data with a power law of the form 
a z oc t a . We made the same analysis for a spiral arm 
mass of M ar mslMdisk = 0.03, an interesting outcome is 
that the value of a is independent of the spiral arm mass. 
Different masses for the spiral arms will just change the 
proportionality constant in the relation a z oc t a . Conse¬ 
quently, a depends only on the pitch angle, and for the 
angles used in our simulation a varies within the range 
0.27 — 0.56 for Sc galaxies. 

Although we have presented the time evolution of a z 
at R = Rmax-, it is possible to made the measure at 
any value of R. Figure Q31 shows the same analysis at 
a different radial position, and this time corresponds to 
the half of the arm length, R = 4.2 kpc. We notice that 
measuring the time evolution of a z at different radius 
give us the same power law behavior. 

With a similar analysis for Sb galaxies, where R max = 
2.5 kpc, we find again that a power law provides a nice 
fit to the data. Figure fl5l shows the er z temporal evolu¬ 
tion for pitch angles of 36° and 45°, and the best fits are 
reached with values for a of 0.32 and 0.37, respectively. 
Smaller angles than those give us plots with more scat¬ 
tered points where a power law fit is not straightforward 
to obtain, i.e., we are not able to establish the value of a 
for angles smaller than 36°. For Sb galaxies we can only 
provide an upper bound for a of ~ 0.37. 


Figure 12. log-log plot of the time evolution of cr z . 


Since a log-log plot is useful to recognize a possible 
power law relationship, Figure [12] shows the log(t) - 
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Figure 13. Temporal evolution of <j z for a Sc galaxy. The panels 
correspond to four different pitch angles, from top to bottom: 50°, 
40°, 30°, and 20°. The black lines show the best fit to the data for 
a relation of the form <j z oc t a . 



Figure 14. Time evolution of a z measured at half arm length for 
a Sc galaxy, 50° pitch angles. The black line shows the best fit to 
the data for a relation of the form cr z oc t a . 



Figure 15. Temporal evolution of cr z for an Sb galaxy. The panels 
correspond with two different pitch angles, 45° (top) and 36° (bot¬ 
tom). The black lines show the best fit to the data for a relation 
of the form a z oc t a . 

Considering that the pitch angles we are employing 
here represent the maximum plausible values of pitch an- 
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gle and spiral arms mass for each type of galaxy, before 
chaos dominates, the values of a presented here would 
represent an upper bound for the contribution of the spi¬ 
ral arms to the vertical velocity dispersion of stars in 
each galaxy type. 

3.1.5. The velocity ellipsoid 

For the the dynamical evolution analysis of the stellar 
disk we applied a classic method derived from the dis¬ 
tribution of the velocity dispersions. The axes of such 
distribution define the known as the “velocity ellipsoid”, 
and this is characterized by the two axes ratios: o z /on 
and o^/on- 

The shape of the velocity ellipsoid seems to show 
a trend in o z Ior with the morphological type 
(jGerssen fc Shapiro Griffinll2012f) . Given the three types 
of galaxies we are using in this work and the parameters 
given in Table [T] that defines them, we are able to mea¬ 
sure the shape of the velocity ellipsoid and compare it 
with the morphological type. In Figure [TG] we plot this 
ratio as measured at the final stage of our simulations. 
The value o z / on decrease with galaxy type and the trend 
keeps with time. 
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Figure 16. Velocity ellipsoid ratio (Tz/ctr as a function of morpho¬ 
logical galaxy type. The dashed line represents a linear regresion 
of the data. 

In Figure [17] we follow the temporal evolution for the 
shape of the velocity ellipsoid. We see that for Sa and Sb 
galaxies o z is always greater than on and the presence of 
spiral arms do not seem to alter this tendency. The initial 
fall in the curves indicates that although o z and on both 
increase with time, on does it in a greater proportion 
than o z , especially for late type galaxies. This heating 
in the i?-direction is similar to that seen in Figure [2] 

The shape of the velocity ellipsoid evolves towards 
lower values and settles down to be nearly constant and 
fluctuates around some equilibrium val ue, a behavior al- 
ready noticed in other numerical work I Sellwoodl[2008( 1. 

Looking at the values of the plots in Figure [TT] we see 
that the velocity dispersions ratio falls to half its initial 
value for Sa and Sb galaxies, with cr 2 always greater than 
on, but the final value of o z /on falls lower than one for 
Sc galaxies. This indicates that the relative increase of 
on for Sc galaxies is greater than that seen on Sa and 


Sb galaxies. This means that the heating effect of the 
spiral arms in Sc galaxies, apart from being more noto¬ 
rious in the z-direction, is also greater in the i?-direction 
compared to the other galaxy types. 



Figure 17. Velocity ellipsoid shape at 5 Gyr time evolution in our 
simulations for the three galaxy types simulated. 


In previous sections we have shown that the vertical 
heating of the stellar disk strongly depends on the pitch 
angle of the spiral pattern. Consequently, the disk thick¬ 
ness in our simulations grows with the value of the pitch 
an gle as seen in Figure [5] a nd equation [TU] As showed 
in iJenkins fe Binnevl (I1990T) . where the ratio o z /on de¬ 
pends on the spiral structure, here we are able to measure 
the shape of the velocity ellipsoid and its dependence on 
the pitch angle. Because the ratio of the velocity disper¬ 
sion reaches a nearly constant value only after a certain 
period of time, under the influence of the spiral pattern 
('Figure fTTl) . we computed the ratio o z /on at the final 
stage in our simulations (5 Gyr) for several pitch angles, 
for the Sc galaxy. In Figure [15] we show how the shape of 
the velocity ellipsoid would be as a function of the pitch 
angle. It decreases first as the pitch angle grows but after 
a 40° angle it fluctuates around some constant value. 

For small pitch angles (< 20°), the radial and vertical 
dispersions are nearly the same and this relation contin¬ 
ues after 5 Gyr, despite the presence of the spiral pattern. 
This is because even when o z and on increase with time, 
the grow rate is small in both directions compared with 
the initial dispersions. On the other hand, Figure [TS] 
shows that for larger pitch angles, not only the incre¬ 
ment in o z is considerable, as it is found in section Rl. 1.11 
but the heating in the i?-direction is remarkably, even 
greater than that in the z-direction. This is expected 
from a structure that rotates in the plane of the galaxy, 
but here we have found that it has a considerable effect 
in both, the radial and vertical directions, and variates 
sensibly with the parameters of the spiral pattern. 
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Figure 18. Velocity ellipsoid shape, measured after a 5Gyr time 
evolution as a function of the pitch angle in an Sc galaxy. 

3.1.6. Transient spiral arms 

All results presented in the previous sections are ob¬ 
tained from simulations that assume long-lived spiral 
arms. However, it is a result from N-body simulations, 
that the s piral arms might rather be transien t features 
in general (ISellwoodll201 It [Kawata et al.ll2011lh 

In this section we present some experiments with tran¬ 
sient spiral arms now, to quantify its effect in the vertical 
structure of the stellar disk and to study possible differ¬ 
ences with our previous results. The adjustability of the 
galaxy models allows us to emulate transient and recur¬ 
rent spiral arms by making them grow and disappear pe¬ 
riodically. The lifetime of the spiral patterns is hard to 
determine in N-body simulations, and even in the cases 
where it can be determined with Fourier modes analysis, 
it is different for every simulation. We take two different 
periods for a simple experiment, one w ith 100 Myr and 
other with 500 Myr ([Grand et al.H201^l . 

The simulations are made with an Sc galaxy model 
with a spiral arms pitch angle of 50°, M arrns /Mdi S k = 
0.05, and two different lifetimes for the transient spiral 
arms. Figure [T9] shows the vertical velocity dispersion 
o z as a function of R for the 5 Gyr evolution time with 
the transient spiral arms. Even with this kind of spiral 
pattern that is not always present in the disk, grows and 
disappear periodically, we can see that <j z keeps increas¬ 
ing with time for both lifetimes used, as it does in our 
previous simulations with longer lasting spiral arms. 

In Figure [20] we compare the final vertical velocity dis¬ 
persion for both models of transient spirals against that 
obtained for non-transient spiral arms. The plots show 
only slight deviations between them, and not significative 
at all radii. 

For a more detailed comparison we pick the radial zone 
at which the increase in vertical velocity dispersion is 
highest and plot the temporal evolution of a z for the 
transient and non-transient spiral arms and show it in 
Figure [7T] Some slight differences are visible at the end 
of the simulation, and as this radial zone is where Aoy 
is maximum, the differences are even smaller for the rest 
of the disk. 

Despite the different nature of the three types of spiral 


2 time(Gyr) 3 



Figure 19. Vertical velocity dispersion as a function of R 
along 5 Gyr evolution time for transient spiral arms with lifetimes 
100 Myr and 500 Myr. 



Figure 20. A comparison of cr z after 5 Gyr for transient spiral 
arms, 100 Myr lifetime, 500 Myr lifetime, and non-transient spiral 
arms. 
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Figure 21. Comparing the time evolution of a z for transient spi¬ 
ral arms, 100 Myr lifetime, 500 Myr lifetime, and non-transient spi¬ 
ral arms. 

arms simulated, the induced stellar dynamical behavior 
is similar. This is because the transient spiral arms, al¬ 
though not always present in the disk, are recurrent and 
form very quickly once the previous spiral pattern dis¬ 
appear as seen in N-body simulations, leaving the stars 
under an almost constant influence. 

3.2. Plausible Origin of the Vertical Stellar Heating 

Based on the set of experiments presented here, we 
conclude that a spiral structure can excite, considerably, 
velocities in the z-direction. Furthermore, we noticed 
that the physical mechanism causing the heating is dif¬ 
ferent from simple resonant excitation. The spiral pat¬ 
tern induces chaotic behavior not linked necessarily to 
resonances, but rather to direct scattering of disk stars, 
which leads to an increase of the velocity dispersion. 

In order to produce evidence to support this gravita¬ 
tional scattering interpretation, we performed the fol¬ 
lowing analysis. For an Sc galaxy with spiral mass 
M arms = 0.05 Mdisk and a pitch angle of 50°, we took a 
sample of 458 particles that at the end of the simulation 
are part of the hot component (i.e. particles that experi¬ 
enced a significant increase in their velocity dispersion). 
By tracing back the initial positions, we reconstruct the 
orbit of each particle. To classify the orbits and dis¬ 
tinguish them as regular or irregular, we im p lemented a 
known spectral method bv lCarpintcro fe Aguilail (ll998ll . 
The method is able also to identify loop, box, and other 
resonant orbits. The results are presented in Figure [22] 

This analysis shows that the orbital sample is com¬ 
pletely dominated by irregular orbits, with no evidence 
of resonances from any kind. In this way we interpret 
the vertical heating as caused by gravitational scattering 
of stars by spiral arms, as already has been pointed out 
in other studies dFuiii et al.ll2011h . 

3.3. The case of the Milky Way 

Probably, the hardest case to study due to the plethora 
of information available, sometimes contradictory, some 
other model or observations dependent, is the Milky Way 
Galaxy. The existence of different types of disks, or dif¬ 
ferent types or vertical structures (the old and young 
thin disks, and the thick disk) has been known but not 


for long. What are the mechanisms that produce the 
heating of stellar orbits into the thick disk?, what is the 
relative importance of each one? is this present in all 
galaxies? the mechanisms that affect radially the disks 
are effective in the vertical component? 

A lot of the finest work on mechanisms to heat dynam¬ 
ically the disk has been done in the plane with steady 
two-dimensional potential models studying fundamen¬ 
tal physical phenomena such as the radial migration. 
Regarding the vertical structure, this is a subject that 
started several decades ago to be of interest in astro¬ 
physics. Only recently, thanks to the calculation power 
of the new generation of supercomputers, we have been 
able to produce more realistic and detailed simulations 
with test particles in steady potentials and with alive 
models with improved resolution, all to study with un¬ 
precedent details the dynamics of the Galactic disk. 

In this section, as a preliminary application to a spe¬ 
cific galaxy, we have constructed a detailed density dis¬ 
tribution, based model of the Milky Way Galaxy that 
includes spiral arms, with the purpose of exploring their 
isolated effect on the heating to the disk, and its par¬ 
ticipation in the process of production of the thick disk. 
We compare our results with o ther recent work d edicated 
to Milky-Way-like potentials ([Faure et al.112014]) . Addi¬ 
tionally, the experiments presented in this section for the 
Milky Way, include one more experiment incorporating 
to the calculations the Galactic bar. We are not pre¬ 
tending this to be an extensive study of the Milky Way 
Galaxy, but rather a modest first approximation and pre¬ 
liminary results of only the vertical heating effects by 
spiral arms and bar on a detailed model of the Galaxy, 
in consequence it is important to mention, we are not in¬ 
cluding all the relevant references to this problem. In an 
ongoing work, we will provide a work fully dedicated to 
the Milky Way Galaxy in this context and the relevant 
references. 

The parameters of the galactic models employed here, 
allow us to model different morphological types and dif¬ 
ferent spiral arms classes by changing the pitch angle, 
the mass or even making them transient. In this sec¬ 
tion, we produce observationally motivated models of the 
Milky Way galaxy and study the response of the stellar 
disk to the spiral arms. The parameters that describe 
the components that compose the axisymmetric back¬ 
ground potential, are constrained by recent estimates of 
the galactic rotation curve. We adopt a pitch angle of 
15.5° for the spiral arms, a mass M arm /Mdi S k = 0.03, 
and an initial time of 0.5 Gyr to increase up to their full 
mass the spiral arms. For the first simulation presented 
for the Milky Way, we do not include a bar in order to 
try to isolate the influence of the spiral pattern in the 
stellar disk (and to compare with other similar work). 

Figure[23](top panel) shows the initial and final vertical 
velocity dispersion after a 5 Gyr evolution for the Milky- 
Way-like model (halo, disk, bulge and spiral arms). As 
seen for the case of early type galaxies (i.e. small pitch 
angles and spiral arm masses), the change in velocity 
dispersion is small. Indeed, although the effect is not 
negligible, in the Milky-Way-like model experiment, with 
a pitch angle of 15.5° and a mass M arm /Mdi S k = 0.03, 
the spiral arms are not capable to heat the stellar disk. In 
Figure [23] (bottom panel) we plot the difference between 
<j z initial and final, and show that the main heating is 
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Figure 22. Orbital classification on the disk particles that experienced a significant increase in their velocity dispersion. The histogram 
shows that the orbital sample is dominated by irregular orbits. 


Table 3 

Parameters for the Milky Way Model 


Parameter 

Value 

Reference 

Axisymmetric Components 

Hq 

8.5 kpc 

1 

00 

220 km s _1 

1 

M b 

1.41 X 10 10 Mq 

1 

M d 

8.56 x 10 10 M 0 

1 

M h 

80.02 x 10 10 M 0 1 

1 

Disk scale-length 

2.5 kpc 

2 

bi 2 

0.3873 kpc 

1 

&2 2 

5.3178 kpc 

1 

b 2 2 

0.2500 kpc 

1 

a3 2 

12 kpc 

1 

Spiral Arms 

locus 

Logarithmic 

3 

arms number 

2 

4,5 

pitch angle 

15.5° 

4 

Msp/Mf) 

3% 


scale-length 

2.5 kpc 

disk based 

Patter speed (Q sp ) 

-20 km s -1 kpc -1 

6 

inner limit 

3.3 kpc 

based on 1LR 

outer limit 

12 kpc 

based on CR 

Triaxial Bar 

major axis 

3.5 kpc 

2,7 

Scale length 

1.7, 0.64, 0.44 kpc 

2 

Axial ratio 

0.64/1.7, 0.44/1.7 


Mass 

1.41 x 10 10 M 0 


Pattern speed (fls) 

50 km s -1 kpc -1 

6 


References. — ( 1 j Allen fe S ant i 11 an 1 9911 (2) |Freudenreichiri998l 
(3) Seigar 8 z James 1998; Seigar et al. 2006. (4) Drimmel et al. 2000. 
(5) Grosb 0 l et al. 2002; Churchwell et al. 2009; Elmegreen 8 z Elmegreen 
2014. (6) Gerhard 2011. (7) Binney et al. 1997; Bissantz Sz Gerhard 
2002. 

Up to 100 kpc halo radius. 

^ hi, a 2 , b 2 , and a .3 are scale lengths. 
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Figure 23. Effect of the Milky-Way-like Spiral Arms in the dy¬ 
namical heating of the disk. Top: Initial and final vertical velocity 
dispersions after 5 Gyr evolution. Bottom: t\(7 z (final-initial). The 
vertical heating is minimum over a small radial region. 

produced in the innermost region of the disk. 

With this experiment we conclude that the vertical 
heating that spiral arms are able to produce in our Milky- 
Way-like Galaxy model is minimum and it is located in a 
small radial region. T his res u lt represen ts a different one 
compared to that of lFaure et al.1 ( 20141 ). that find some 
important effect on the vertical heating due to the spiral 
arms. To try to account for the origin of this discrepancy, 
further studies will be performed in a future work. 

In a second experiment, we have also included a known 
mechanism that exerts strong secular radial and vertical 
evolution on galactic disks, the central bar. Although 
the present work is advocated to the study of the effect 
of the spiral arms on stellar disks vertical heating, we 
have produced a careful experiment for the Milky Way 
Galaxy to explore to what extent the central bar can 
contribute to the vertical heating. In an ongoing paper 
we study in detail the vertical heating produced by the 
full galactic model of the Milky Way, here we present 
some preliminary results of this study. 

We ran a simulation with the parameters representa¬ 
tive of the Milky Way Galaxy, given on Table [3j Fig¬ 
ure [24] shows the initial and final vertical velocity dis¬ 
persion <j z across the stellar disk. Although the two 
nonaxisymmetric structures are present in the simula¬ 
tion, by comparing Figures E51 and UHl there is no doubt 
that the heating seen in the second one is produced by 
the galactic bar, and that the only affected region seems 
to be the central part (sim ilar results were obtained by 
(ISaha. Tseng fc Taamll2010t )l. 

The effect of the galactic bar on the vertical velocity 
dispersion is considerably, with an increase in a z up to 
16kms _1 , mostly in the inner 6kpc of the disk. 

4. CONCLUSIONS 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

R(kpc) 


Figure 24. Effect of the Milky Way Spiral Arms + Bar model on 
the vertical disk Heating. Top: Initial and final vertical velocity 
dispersions after 5 Gyr evolution. Bottom: A<r z (final-initial). The 
vertical heating is appreciable, with an increment up to 16kms —1 
in the inner region of the disk. 

With the use of detailed models for spiral galaxies, we 
produce an extensive study of disks vertical heating. The 
models include an axisymmetric component (bulge, mas¬ 
sive halo and disk) plus a density based three dimen¬ 
sional potential for the spiral arms, orbitally tested for 
self-consistency. 

The main outcome of this study is that the spiral 
structure is capable to excite moderate to high disper¬ 
sion velocities in the z-direction inducing a relation be¬ 
tween stellar age and u z in all spirals except the earliest 
types. Although rather small, this mechanism works for 
the Milky Way Galaxy, assuming the parameters for the 
spiral arms known to this day. Consequently, by isolat¬ 
ing the effect of the spiral arms on the vertical velocity 
dispersion of the stellar disk, we can conclude that spi¬ 
ral arms have the capability to contribute to the heating 
mechanism that gives raise to thick disks in spiral galax¬ 
ies, from intermediate to late morphological types. 

Therefore, the thickness of the stellar disk driven by 
the spiral arms go from negligible to very important, de¬ 
pending on the characteristics of the spiral pattern such 
as mass, pitch angle and angular speed along with the 
morphological galaxy type. By covering a whole set of 
values for this parameters in the test particle simulations 
we conclude that: 

• The relative increase in vertical velocity disper¬ 
sion is a function of the morphological type, be¬ 
ing smaller for early type galaxies (Sa type in the 
Hubble scheme) and larger from intermediate to 
late type galaxies. Although for the sake of clarity, 
we present our results by separating the models in 
approximately average Sa, Sb and Sc galaxy types, 
we are aware that strong gaps are in between the 
different models in the initial scale height. The 
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results on the disk heating driven by spiral arms 
presented here, are therefore more general, i.e. the 
observed heating is considerably more significant in 
thinner disks, which, as a consequence, has impli¬ 
cations on galaxy types, in this case, particularly 
on later types. 

Two distinctive characteristics of the spiral arms 
are the mass and the pitch angle. By varying this 
parameters one at the time, within the permitted 
values, we found that they play equally important 
roles in the vertical heating of the stellar disk. As 
expected, massive arms and/or large pitch angles 
are responsible for the most prominent disk thick¬ 
ness found in our simulations. 

The other dynamical parameter characteristic of 
the spiral arms, is the pattern speed f2. When vary¬ 
ing this value, we found that the vertical heating 
induced by the spiral pattern is greater for slow 
rotating arms and decrease for those that rotate 
faster. Smaller values of Q allow the spiral arms to 
interact more effectively with the stars, and heat 
more efficiently the stellar disk, independently of 
galaxy type. 

By measuring the time evolution of a z in our simu¬ 
lations, we found that the vertical heating induced 
by the spiral arms follows a power law a z oc f“, as 
seen from observations in the solar neighbourhood 
and other numerical simulations. An interesting 
result is that the value of a depends only on the 
pitch angle and not on the mass of the arms. 

The ratio ct z /ctr, that defines the shape of the ve¬ 
locity ellipsoid, shows a clear trend with the mor¬ 
phological type that remains during the temporal 
evolution. The values of this ratio across the en¬ 
tire simulations show that the heating in the In¬ 
direction is remarkable, always greater that in the 
z-direction for the three galaxy types. This is ex¬ 
pected from a structure that rotates in the plane 
of the galaxy, but here we have found that it has a 
considerably effect in both, the radial and vertical 
directions. Consequently, the spiral arms are key 
to determine the shape of the velocity ellipsoid. 

Although in this work we cover the general prop¬ 
erties of normal spiral galaxies, we took advantage 
of the adjustability of our model to represent the 
gravitational potential of the Milky Way Galaxy. 
Analysing the change in velocity dispersion induced 
by the nonaxisymmetric structures of the Galaxy 
we found that: the galactic spiral arms are not ca¬ 
pable to induce an important thickness in the stel¬ 
lar disk, the increase in vertical velocity dispersion 
is small, therefore, from this study we conclude 
that the spiral arms play no role in producing a 
thick disk. If we add the galactic bar on the other 
hand, the vertical velocity dispersion increase con¬ 
siderably, mostly within the region covered by the 
bar. This means that for the Milky Way, the bar 
is an important heating mechanism that should be 
considered in calculations, but mostly in the inner 
region of the disk. 
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